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We discuss one of the most interesting phenomena exhibited by baby skyrmions - 
breaking of rotational symmetry. The topics we will deal with here include the ap- 
pearance of rotational symmetry breaking in the static solutions of baby Skyrme 
models, both in flat as well as in curved spaces, the zero-temperature crystalline 
structure of baby skyrmions, and finally, the appearance of spontaneous breaking 
of rotational symmetry in rotating baby skyrmions. 

& . 

. 1.1. Breaking of Rotational Symmetry in Baby Skyrme Models 



The Skyrme model 1,2 is an SU(2)-valucd nonlinear theory for pions in (3+1) di- 
mensions with topological soliton solutions called skyrmions. Apart from a kinetic 
term, the Lagrangian of the model contains a 'Skyrme' term which is of the fourth 
order in derivatives, and is used to introduce scale to the model. 3 The existence 
of stable solutions in the Skyrme model is a consequence of the nontrivial topol- 
ogy of the mapping M of the physical space into the field space at a given time, 
O \ M : S 3 — > SU(2) = S 3 , where the physical space R 3 is compactified to S 3 

by requiring the spatial infinity to be equivalent in each direction. The topology 
which stems from this one-point compactification allows the classification of maps 
into equivalence classes, each of which has a unique conserved quantity called the 
^ . topological charge. 

| The Skyrme model has an analogue in (2+1) dimensions known as the baby 

Skyrme model, also admitting stable field configurations of a solitonic nature. 4 Due 
to its lower dimension, the baby Skyrme model serves as a simplification of the 
original model, but nonetheless it has a physical significance in its own right, hav- 
ing several applications in condensed-matter physics, 5 specifically in ferromagnetic 
quantum Hall systems. 6-9 There, baby skyrmions describe the excitations relative 
to ferromagnetic quantum Hall states, in terms of a gradient expansion in the spin 
density, a field with properties analogous to the pion field in the 3D case. 10 

The target manifold in the baby model is described by a three-dimensional vector 
4> = (0i, 02, 03 ) with the constraint </>• <fi = 1. In analogy with the (3 + 1)D case, the 
domain of this model R 2 is compactified to S 2 , yielding the topology required for 

1 to appear in: G. Brown and M. Rho, Eds., Multifaceted Skyrmions, (World Scientific, Singapore, 
2009). 

2 A version of this manuscript with higher-resolution figures is available at 
www . tau. ac . il/~itayhe/SkReview/SkReview. rar 
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the classification of its field configurations into classes with conserved topological 
charges. The Lagrangian density of the baby Skyrme model is given by: 

C = ■ 5*0 - y [{d^ ■ d^f - (drf ■ d v ct>) ■ (d^cf, ■ <?»] - Ufa) , 

(1.1) 

and consists of a kinetic term, a Skyrme term and a potential term. 

While in (3+1) dimensions the latter term is optional, 11 its presence in the 
(2+l)D model is necessary for the stability of the solutions. However, aside from 
the requirement that the potential vanishes at infinity for a given vacuum field value 
(normally taken to be = (0, 0, 1)), its exact form is arbitrary and gives rise to 
a rich family of possible baby-Skyrme models, several of which have been studied 
in detail in the literature. The simplest potential is the 'holomorphic' model with 
U(4>3) = /i 2 (l — (fe) 4 - 12-14 It is known to have a stable solution only in the charge- 
one sector (the name refers to the fact that the stable solution has an analytic 
form in terms of holomorphic functions). The model with the potential t/(</>3) = 
/i 2 (l — 03) (commonly referred to as the 'old' model) has also been extensively 
studied. This potential gives rise to very structured non-rotationally-symmetric 
multi-skyrmions. 4,15 Another model with U{4>-z) = /j, 2 (1 — 4>\) produces ring-like 
multi-skyrmions. 16 Other double- vacuum potentials which give rise to other types 
of solutions have also been studied. 17 

Clearly, the form of the potential term has a decisive effect on the properties 
of the minimal energy configurations of the model. It is then worthwhile to see 
how the multisolitons of the baby Skyrme model look like for the one-parametric 
family of potentials U = fJ. 2 (l — 4>3) s which generalizes the 'old' model (s = 1) and 
the holomorphic model (s = 4). 18 As it turns out, the value of the parameter s 
has dramatic effects on the static solutions of the model, both quantitatively and 
qualitatively, in the sense that it can be viewed as a 'control' parameter responsible 
for the repulsion or attraction between skyrmions, which in turn determines whether 
or not the minimal-energy configuration breaks rotational symmetry. 

The Lagrangian density is now: 

c = ■ d»<p - y ({drf ■ d^f - (d,4> ■ d v <t>) ■ [d»<t> ■ d v 4>)) - m 2 (i - 3 ) s , 

(1.2) 

and contains three free parameters, namely K, [i and s. Since either k or /i may be 
scaled away, the parameter space of this model is in fact only two dimensional. Our 
main goal here is to study the effects of these parameters on the static solutions of 
the model within each topological sector. 

The multi-skyrmions of our model are those field configurations which minimize 
the static energy functional within each topological sector. In polar coordinates the 
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energy functional is given by 

E = Jrdrd9 Q(9 r • d r <f> + ^d 8 4> ■ decP) + ^Al^t + ^(i _ ^ . 

(1.3) 

The Euler-Lagrange equations derived from the energy functional (|1.3p are non- 
linear PDE's, so in most cases one must resort to numerical techniques in order to 
solve them. In our approach, the minimal energy configuration of a baby skyrmion 
of charge B and a given set of values /i, «, s is found by a full-field relaxation method, 
which we describe in more detail in the Appendix. 

1.1.1. Results 

Applying the minimization procedure, we obtain the static solutions of the model for 
1 < B < 5. Since the parameter space of the model is effectively two dimensional 
(as discussed earlier), without loss of generality we fix the potential strength at 
fi 2 = 0.1 throughout, and the s-k parameter space is scanned in the region < s < 4, 
0.01 < k 2 < 1. 

1.1.1.1. Charge-one skyrmions 

In the charge-one sector, the solutions for every value of s and n are stable 
rotationally-symmetric configurations. Figure [TTTk shows the obtained profile func- 
tions of the B = 1 solution for different values of s with k fixed at k 2 = 0.25. 
Interestingly, the skyrmion energy as a function of s is not monotonic, but acquires 
a minimum at s ~ 2.2, as is shown in Fig. 11.21 




123456 12345S 



(a) B=l (b) B=2 



Fig. 1.1. Profile functions of the B = 1 (left) and B = 2 (right) skyrmions for s = 0.5 (solid), 
s = 1 (dotted) and s = 2 (dot-dashed). Here re is fixed at re 2 = 0.25. 
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Fig. 1.2. Total energies (divided by AnB) of the charge-one (♦) charge-two (I) and charge-three 
(0) skyrmions as a function of the parameter s for k 2 = 0.05. Each of the energy graphs attains 
a minimal value at some s. At s 2 the energy-per-topological-charge of the charge- two and 
charge-three solutions reaches the charge-one energy (from below), and stable solutions are no 
longer observed. 



1.1.1.2. Charge-two skyrmions 

Stable solutions also exist in the B = 2 sector, but only for s < 2. They are 
rotationally-symmetric and ring-like, corresponding to two charge-one skyrmions 
on top of each other. Figure ITTTb shows the profile functions of the stable solutions 
for different values of s and k 2 = 0.25. 

As in the B = 1 case, the energy of the charge-two skyrmion as a function of 
s is non-monotonic and has a minimum around s = 1.3. As shown in Fig. 11.21 at 
s ss 2 the energy of the ring-like configuration reaches the value of twice the energy 
of the charge-one skyrmion and stable configurations cease to exist. At this point, 
the skyrmion breaks apart into its constituent charge-one skyrmions, which in turn 
start drifting away from each other, thus breaking the rotational symmetry of the 
solution. Contour plots of the energy distribution of the B = 2 skyrmion arc shown 
in Fig. 11.31 for k 2 = 1 and for two s values. While for s = 1.5 a ring-like stable 
configuration exists (Fig. 11.3b ). for s = 2.6 the skyrmion breaks apart. The latter 
case is shown in Fig. 11.3b which is a "snapshot" taken while the distance between 
the individual skyrmions kept growing. 

These results are in accord with corresponding results from previously known 
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studies of both the 'old' (s = 1) model in which ring- like configurations have been 
observed, 4,15 and the holomorphic model for which no stable solutions have been 
found. 12,13 




(a) s = 1.5 (b) s = 2.6 



Fig. 1.3. Contour plots of the energy distributions (ranging from violet - low density to red — 
high density) of the B = 2 skyrmion for re 2 = 1. In the s < 2 regime, ring-like rotationally- 
symmetric configurations exist, corresponding to two charge-one skyrmions on top of each other 
(left), whereas for s > 2, the charge-two skyrmion splits into two one-charge skyrmions drifting 
infinitely apart (right). 

Rotationally-symmetric charge- two locally stable solutions may also be observed 
in the large s regime, including the 'holomorphic' s = 4 case, in which case the 
global minimum in this regime corresponds to two infinitely separated charge-one 
skyrmions. The total energy of the rotationally symmetric solutions is larger than 
twice the energy of a charge-one skyrmion, indicating that the split skyrmion is an 
energetically more favorable configuration. We discuss this issue in more detail in 
the section [1. 2. 2 1 

1.1.1.3. Charge-three and higher-charge skyrmions 

As with the B = 2 skyrmion, the existence of stable charge-three skyrmions was also 
found to be s dependent. For any tested value of k in the range 0.01 < k 2 < 1, we 
have found that above s as 2, no stable charge-three solutions exist; in this region 
the skyrmion breaks apart into individual skyrmions drifting further and further 
away from each other. 

In the s < 2 region, where stable solutions exist, the energy distribution of the 
charge-three skyrmion turns out to be highly dependent on both s and k. Keeping 
k fixed and varying s, we find that in the small s regime, ring-like rotationally- 
symmetric configurations exist. Increasing the value of s results in stable minimal 
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energy configurations with only Z(2) symmetry, corresponding to three partially- 
overlapping charge-one skyrmions in a row. as shown in Figs 11.4b and 11.4b . The 
energy of the charge-three skyrmion also has a minimum in s, at around s = 1.5 (as 
shown in Fig. II. 2|l . At s 2 the energy of the skyrmion becomes larger than three 
times the energy of a charge-one skyrmion and stable configurations arc no longer 
obtainable. This is illustrated in Fig. 11.41 which shows contour plots of the energy 
distribution of the B = 3 skyrmion for different values of s and fixed k. For s = 0.5 
(Fig. 11.4b ). the solution is rotationally symmetric and for s = 0.75 and s = 1 (Figs 
11.4b and 11.4b respectively) the rotational symmetry of the solution is broken and 
only Z(2) symmetry remains. At s = 3, no stable solution exists. The latter case 
is shown in Fig. I1.4b l which is a "snapshot" taken while the distance between the 
individual skyrmions kept growing. 

The dependence of the skyrmion solutions on the value of n with fixed s show 
the following behavior: While for small k the minimal energy configurations are 
rotationally-symmetric, increasing its value results in an increasingly larger rota- 
tional symmetry breaking. This is illustrated in Fig. 11.51 




case, the minimal energy configuration is rotationally symmetric, corresponding the three one- 
skyrmions on top of each other. For s = 0.75 and s = 1 the solutions exhibit only 2(2) symmetry, 
corresponding to partially-overlapping one-skyrmions. For s = 3 no stable solution exists and the 
individual skyrmions arc drifting apart. 
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(a) re 3 = 0.01 (b) k 2 = 0.1 (c) K 2 = 0.25 (d) f,: 2 = 1 



Fig. 1.5. Energy densities and corresponding contour plots (ranging from violet - low density to 
red - high density) of the B = 3 skyrmion for fixed s (s = 0.5) and varying re. At low re, the 
minimal energy configuration is rotationally symmetric. As re is increased, breaking of rotational 
symmetry appears, and only Z(2) symmetry remains. 



The B = 4 and B = 5 skyrmion solutions behave similarly to the B = 3 
solutions. This is illustrated in Fig. 11.61 which shows the stable solutions that have 
been obtained in the s = 0.9 case and the splitting of these skyrmions into their 
constituents in the s — 4 case. 



1.2. The Lattice Structure of Baby Skyrmions 

The Skyrme model 1 may also be used to describe systems of a few nucleons, and has 
also been applied to nuclear and quark matter. 19-21 One of the most complicated 
aspects of the physics of hadrons is the behavior of the phase diagram of hadronic 
matter at finite density at low or even zero temperature. Particularly, the properties 
of zero-temperature skyrmions on a lattice are interesting, since the behavior of 
nuclear matter at high densities is now a focus of considerable interest. Within the 
standard zero-temperature Skyrme model description, a crystal of nucleons turns 
into a crystal of half nucleons at finite density. 22-26 

To study skyrmion crystals one imposes periodic boundary conditions on the 
Skyrme field and works within a unit cell. 11 The first attempted construction of a 
crystal was by Klebanov, 22 using a simple cubic lattice of skyrmions whose symme- 
tries maximize the attraction between nearest neighbors. Other symmetries were 
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(a) B=4, s=0.9 



(b) B=5, s=0.9 




(c) B=4, s=4 



(d) B=5, s=4 



Fig. 1.6. Contour plots of the energy distributions (ranging from violet - low density to red — 
high density) of the B = 4 and B = 5 skyrmions for s = 0.9 and s = 4 (k 2 = 0.1). In the lower 
s region stable solutions exist; the upper figures show a B = 4 skyrmion in a bound state of two 
charge- two skyrmions (left), and a B = 5 skyrmion in a two-one- two configuration. For values of 
s higher than 2, the multi-skyrmions split into individual one-skyrmions constantly drifting apart 
(lower figures). 



proposed which lead to crystals with slightly lower, but not minimal energy. 23,2 It 
is now understood that it is best to arrange the skyrmions initially as a face-centered 
cubic lattice, with their orientations chosen symmetrically to give maximal attrac- 
tion between all nearest neighbors. 25 ' 26 

The baby Skyrme model too has been studied under various lattice settings 27-31 
and in fact, it is known that the baby skyrmions also split into half-skyrmions when 
placed inside a rectangular lattice. 29 However, as we shall see, the rectangular 
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periodic boundary conditions do not yield the true minimal energy configurations 
over all possible lattice types. 32 This fact is particularly interesting both because of 
its relevance to quantum Hall systems in two-dimensions, and also because it may 
be used to conjecture the crystalline structure of nucleons in three-dimensions. 

In two dimensions there are five lattice types, as given by the crystallographic 
restriction theorem. 33 In in all of them the fundamental unit cell is a certain type of 
a parallelogram. To find the crystalline structure of the baby skyrmions, we place 
them inside different parallelograms with periodic boundary conditions and find the 
minimal energy configurations over all parallelograms of fixed area (thus keeping the 
skyrmion density fixed). As we show later, there is a certain type of parallelogram, 
namely the hexagonal, which yields the minimal energy configuration. In particular, 
its energy is lower than the known 'square-cell' configurations in which the skyrmion 
splits into half-skyrmions. As will be pointed out later, the hexagonal structure 
revealed here is not unique to the present model, but also arises in other solitonic 
models, such as Ginzburg-Landau vortices, 34 quantum Hall systems, 6,7 and even in 
the context of 3D skyrmions. 35 The reason for this will also be discussed later. 

In what follows we review the setup of our numerical computations, introducing a 
systematic approach for the identification of the minimal energy crystalline structure 
of baby skyrmions. In section fl. 2. 21 we present the main results of our study and in 
section [T.2.31 a somewhat more analytical analysis of the problem is presented. 

1.2.1. Baby skyrmions inside a parallelogram 

We find the static solutions of the model by minimizing the static energy functional: 

E= ^J^dxdy((d x <j>) 2 + (d y <{>) 2 +K 2 (d x <f>x cV£) 2 + 2/i 2 (l-0 3 )) , (1.4) 

within each topological sector. In this example, we use the 'old' model potential 
term. In our setup, the integration is over parallelograms, denoted here by A: 

A = {ai(L, 0) + a 2 (sisin7, SLCOS7) : < ai,a^ < 1} . (1.5) 

Here L is the length of one side of the parallelogram, sL with < s < 1 is the 
length of its other side and 0<7<7r/2is the angle between the l sL' side and the 
vertical to the side (as sketched in Fig. II. 7|) . 

Each parallelogram is thus specified by a set {L, s, 7} and the skyrmion density 
inside a parallelogram is p = Bj (sL 2 cos 7), where B is the topological charge of the 
skyrmion. The periodic boundary conditions are taken into account by identifying 
each of the two opposite sides of a parallelogram as equivalent: 

4>{x) = 4>(x + ri\(L, 0) + n 2 (sl/sin7, sLcos7)) , (1.6) 

with rei,7i2 € Z. We are interested in static finite-energy solutions, which in the 
language of differential geometry are T 2 ► S2 maps. These are partitioned into 
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Fig. 1.7. The parameterization of the fundamental unit-cell parallelogram A (in black) and the 
two-torus T2 into which it is mapped (in gray). 

homotopy sectors parameterized by an invariant integral topological charge B, the 
degree of the map, given by: 

B = i / dxdy (cp • {d x (t> x dyej})) . (1.7) 

The static energy E can be shown to satisfy 

£> 471-5, (1.8) 

with equality possible only in the 'pure' 0(3) case (i.e., when both the Skyrme and 
the potential terms are absent). 29 We note that while in the baby Skyrme model 
on IR 2 with fixed boundary conditions the potential term is necessary to prevent 
the solitons from expanding indefinitely, in our setup it is not required, due to the 
periodic boundary conditions. 29 We study the model both with and without the 
potential term. 

The problem in question can be simplified by a linear mapping of the paral- 
lelograms A into the unit-area two-torus Tg. In the new coordinates, the energy 
functional becomes 

E = [ dxdy (s 2 ^) 2 - 2ssm 1 d x (f>d y (f>+ (d y cf>) 2 ) 

2s cos 7 Jt 2 

+ tw f dxd y x d v<i>? + — [ dxd v (i - <fe) • (i.9) 

We note that the dependence of the energy on the Skyrme parameters k and \i and 
the skyrmion density p is only through k 2 p and (J? /p. 
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In order to find the minimal energy configuration of skyrmions over all paral- 
lelograms with fixed area (equivalcntly, a specified p), we scan the parallelogram 
parameter space {s, 7} and find the parallelogram for which the resultant energy is 
minimal over the parameter space. An alternative approach to this problem, which 
is of a more analytical nature, may also been implemented. We discuss it in detail 
in section [1. 2. 31 

1.2.2. Results 

In what follows, we present the minimal energy static skyrmion configurations over 
all parallelograms, for various settings: The 'pure' 0(3) case, in which both n, the 
Skyrme parameter, and p, the potential coupling, are set to zero, the Skyrme case 
for which only p = 0, and the general case for which neither the Skyrme term nor 
the potential term vanish. 

In each of these settings, we scanned the parameter space of parallelograms, 
while the skyrmion density p was held fixed, yielding for each set of {s,7} a mini- 
mal energy configuration. The choice as to how many skyrmions are to be placed 
inside the unit cells was made after some preliminary testing in which skyrmions of 
other charges (up to B = 8) were also examined. The odd-charge minimal-energy 
configurations turn out to have substantially higher energies than even-charge ones, 
where among the latter, the charge-two skyrmion is found to be the most funda- 
mental, as it is observed that the charge-two configuration is a 'building-block' of 
the higher-charge configurations. This is illustrated in Fig. ll.Sl in which the typical 
behavior of the multi-skyrmion energies as a function of topological charge is shown. 



1.2.2.1. The pure 0(3) case (k = p = 0) 

The pure 0(3) case corresponds to setting both k and p to zero. In this case, 
analytic solutions in terms of Weierstrass elliptic functions may be found 29-31 and 
the minimal energy configurations, in all parallelogram settings, saturate the energy 
bound in l|1.8j) giving E = AttB. Thus, comparison of our numerical results with 
the analytic solutions serves as a useful check on the precision of our numerical 
procedure. The agreement is to 6 significant digits. Contour plots of the charge 
densities for different parallelogram settings for the charge- two skyrmions are shown 
in Fig. 11.91 all of them of equal energy E/8ir = 1. 

1.2.2.2. The Skyrme case ( K ^0,p = 0) 

As pointed out earlier, for p = the dependence of the energy functional on the 
Skyrme parameter n is only through K 2 p, so without loss of generality we vary p 
and fix k 2 — 0.03 throughout (this particular choice for k was made for numerical 
convenience). Minimization of the energy functional (|1.9p over all parallelograms 
yields the following. For any fixed density p, the minimal energy is obtained for 
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Fig. 1.8. Energy per charge of the multi-skyrmion configurations as a function of topological 
charge. The horizontal dashed line was added to guide the eye. (Here, k 2 = 0.03, fj, = 0, p = 1, 
s = 1 and 7 = 7r/6.) 



s = 1 and 7 = 7r/6. This set of values corresponds to the 'hexagonal' or 'equilateral 
triangular' lattice. In this configuration, any three adjacent zero-energy loci (or 
'holes') are the vertices of equilateral triangles, and eight distinct high-density peaks 
are observed (as shown in Fig. II. 10b ). This configuration can thus be interpreted 
as the splitting of the two-skyrmion into eight quarter-skyrmions. This result is 
independent of the skyrmion density p. 

In particular, the well-studied square-cell minimal energy configuration (Fig. 
11.10b ). in which the two-skyrmion splits into four half-skyrmions, has higher energy 
than the hexagonal case. Figure [TTTOl shows the total energies (divided by 87r) and 
the corresponding contour plots of charge densities of the hexagonal, square and 
other configurations (for comparison), all of them with p = 2. 

The total energy of the skyrmions in the hexagonal setting turns out to be 
linearly proportional to the density of the skyrmions, reflecting the scale invariancc 
of the model (Fig. II. lip . In particular, the global minimum of E = 4ttB = 8tt is 
reached when p — ► 0. This is expected since setting p = is equivalent to setting 
the Skyrme parameter k to zero, in which case the model is effectively pure 0(3) 
and inequality (|1.8|) is saturated. 
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(a) s — 1 and 7 — 







(c) s = 0.7 and 7 = 




(b) s = 0.9 and 7 =tt/16 




(d) s = 1 and 7 = tt/6 



Fig. 1.9. Charge-two skyrmions in the pure 0(3) case: Contour plots of the charge densities 
ranging from violet (low density) to red (high density) for various parallelogram settings, all of 
which saturate the energy bound E = 4irB = 87r. 



1.2.2.3. The general case (n ^ 0, p =/= 0) 

The hexagonal setting turns out to be the energetically favorable also in the general 
case. Moreover, since in this case the skyrmion has a definite size (as is demonstrated 
by the p dependence in the energy functional) , the skyrmion structure is different at 
low and at high densities and a phase transition occurs at a certain critical density. 
While at low densities the individual skyrmions are isolated from each other, at 
high densities they fuse together, forming the quarter-skyrmion crystal, as in the 
Skyrme case reported above. As the density p decreases, or equivalently the value 
of p increases, the size of the skyrmions becomes small compared to the cell size. 
The exact shape of the lattice loses its effects and the differences in energy among 
the various lattice types become very small. This is illustrated in Fig. 11.111 

Due to the finite size of the skyrmion, there is an optimal density for which the 
energy is minimal among all densities. Figure IT. 121 shows the contour plots of the 
charge density of the charge-two skyrmion for several densities with k 2 = 0.03 and 
p 2 = 0.1. The energy of the skyrmion is minimal for p f=a 0.14 (Fig. Il.ll|) . 
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(c) * - 0.51, 7 5= 0. and (d) s 
E/Stt = 1.587 



= 1, 7 = 7r/4, and 
B/Stt = 1.454 



Fig. 1.10. Charge-two skyrmions in the Skyrme case with re 2 = 0.03 and p = 2: Contour plots of 
the charge densities for the hexagonal, square and other fundamental cells ranging from violet (low 
density) to red (high density). As the captions of the individual subfigurcs indicate, the hexagonal 
configuration has the lowest energy. 



1.2.3. Semi- analytical approach 

The energy functional (|1.9[) depends both on the Skyrme field <fr and on the paral- 
lelogram parameters 7 and s. Formally, the minimal energy configuration over all 
parallelograms may be obtained by functional differentiation with respect to <fi and 
regular differentiation with respect to 7 and s. However, since the resulting equa- 
tions are very complicated, a direct numerical solution is quite hard. Nonetheless, 
some analytical progress may be achieved in the following way. As a first step, we 
differentiate the energy functional (|1.9|) only with respect to 7 and s: 



dE 



1 



(?7 2 s cos 2 7 
dE _ 1 
ds 2s 2 cos 7 



(smj(S yy + s 2 £ xx )-2s£ xy 



0, 



{£ m - s z £ 



2 exx \ 



0, 



(1.10) 
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Fig. 1.11. Total energy E (divided by 8n) of the charge-two skyrmion in the hexagonal lattice (I 
- Skyrme case and ♦ - general case) and in the square lattice ( ■ - Skyrme case and ♦ - general 
case) as function of the skyrmion density (in the Skyrme case, re 2 = 0.03 and in the general case 
k 2 = 0.03 and p? = 0.1). Note the existence of an optimal density (at p R3 0.14) in the general 
case, for which the energy attains a global minimum. Figure (b) is an enlargement of the lower 
left corner of figure (a). 



where £ l ° = dxdy(di<p ■ dj(p) and i,j £ {x, y}. Solving these equations for 7 and 
s yields 




(1.11) 



£xy 

sin 7 = z . 

Substituting these expressions into the energy functional (|1.9[) . we arrive at a 're- 
duced' functional 




(1.12) 



where £ s k = J T ^ dxdy {d x <fi x d y 4>) 2 is the Skyrme energy and £ pot = 
dxdy (1 — 03) is the potential energy. Now that both 7 and s are eliminated 
from the resultant expression, and the conditions for their optimization are built 
into the functional, the numerical minimization is carried out. We note here, how- 
ever, that the procedure presented above should be treated with caution. This is 
since Eqs. are only cxtremum conditions, and may correspond to a maximum 

or saddle-point. Hence, it is important to confirm any results obtained using this 
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(a) p= 1. E/Sit = 1.220 




(b) p = 0.5. E/8tt = 1.128 (c) p = 0.25. E/Stv = 1 ,087 




(d) p = 0.14. £7/8tt = 1.08 (e) p = 0.0625. £/8tt = 1.082 (f) p = 0.03125, E/8tt = 1.084 



Fig. 1.12. Charge-two skyrmions in the general case with pi 2 = 0.1 and ft 2 = 0.03: Contour plots 
of the charge densities of the minimal-energy configurations in the hexagonal setting for different 
densities. Here, the energetically most favorable density is p = 0.14. The plot colors range from 
violet (low density) to red (high density). 



method by comparing them with corresponding results obtained from the method 
described in the previous section. 

It is therefore reassuring that numerical minimization of the reduced functional 
(|1.12p gives sin 7 = 0.498 (7 « n/6) and s = 1 (both for the Skyrme case and the 
general case), confirming the results presented in the previous section. 

In the general (/x ^ 0) case, the energy functional (|1.12[) may be further differen- 
tiated with respect to the skyrmion density p to obtain the optimal density for which 
the skyrmion energy is minimal. Differentiating with respect to p, and substituting 
the obtained expression into the energy functional, results in the functional 

E = \j£ xx £yy - {E x vf + Kp^2£ sV E pot . (1.13) 

Numerical minimization of the above expression for k 2 — 0.03 and various p values 
(0.1 < p 2 < 10) yields the hexagonal setting as in the Skyrme case. In particular, 
for p 2 — 0.1 the optimal density turns out to be p ~ 0.14, in accord with results 
presented in Sec. 11.2.2.31 

1.2.4. Further remarks 

As pointed out earlier, the special role of the hexagonal lattice revealed here is 
not unique to the baby Skyrme model, but in fact arises in other solitonic models. 
In the context of Skyrme models, the existence of a hexagonal two dimensional 
structure of 3D skyrmions has also been found by Battye and Sutcliffe, 35 where it 
has already been noted that energetically, the optimum infinite planar structure of 
3D skyrmions is the hexagonal lattice, which resembles the structure of a graphite 
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sheet, the most stable form of carbon thcrmodynamically. 11 Other examples in 
which the hexagonal structure is revealed are Ginzburg-Landau vortices which are 
known to have lower energy in a hexagonal configuration than on a square lattice. 34 
Thus, it should not come as a surprise that the hexagonal structure is found to be 
the most favorable in the baby Skyrme model. 

As briefly noted in the opening paragraphs of this section, a certain type of 
baby skyrmions also arise in quantum Hall systems as low-energy excitations of the 
ground state near ferromagnetic filling factors (notably 1 and 1/3). 6 It has been 
pointed out that this state contains a finite density of skyrmions, 36 and in fact the 
hexagonal configuration has been suggested as a candidate for their lattice struc- 
ture. 7 Our results may therefore serve as a supporting evidence in that direction, 
although a more detailed analysis is in order. 

Our results also raise some interesting questions. The first is how the dynamical 
properties of baby skyrmions on the hexagonal lattice differ from their behavior in 
the usual rectangular lattice. Another question has to do with their behavior in 
non-zero-temperature. 

One may also wonder whether and how these results can be generalized to the 
3D case. Is the face-centered cubic lattice indeed the minimal energy crystalline 
structure of 3D skyrmions among all parallelepiped lattices? If not, what is the 
minimal energy structure, and how do these results depend on the presence of a 
mass term? These questions await a systematic study. 



1.3. Baby Skyrmions on the Two-Sphere 

Although skyrmions were originally introduced to describe baryons in three spatial 
dimensions, 1 they have been shown to exist for a very wide class of geometries, 37 
specifically cylinders, two-spheres and three-spheres. 38-42 

Here, we consider a baby Skyrme model on the two-sphere 0. We compute the 
full-field minimal energy solutions of the model up to charge 14 and show that they 
exhibit complex multi-skyrmion solutions closely related to the skyrmion solutions 
of the 3D model with the same topological charge. To obtain the minimum energy 
configurations, we apply two completely different methods. One is the full-field 
relaxation method, with which exact numerical solutions of the model are obtained, 
and the other is a rational map approximation scheme, which as we show yields 
very good approximate solutions. 

In an exact analogy to the 3D Skyrme model, the charge-one skyrmion has a 
spherical energy distribution, the charge-two skyrmion is toroidal, and skyrmions 
with higher charge all have point symmetries which are subgroups of 0(3). As we 
shall see, it is not a coincidence that the symmetries of these solutions are the same 
as those of the 3D skyrmions. 



a This type of model has been studied before, 39,40 although only rotationally-symmetric configu- 
rations have been considered. 
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1.3.1. The baby Skyrme model on the two-sphere 

The model in question is a baby Skyrme model in which both the domain and target 
are two-spheres. The Lagrangian density here is simply 

C = ■ d»ct> + y ■ d»(t>) 2 - (d^ ■ d u <t>){d»<t> ■ 0»] , (1.14) 

with metric ds 2 = dt 2 — d6 2 — sin 2 8 dip 2 , where 8 is the polar angle £ [0, ir] and 
(p is the azimuthal angle £ [0, 2tt). The Lagrangian of this model is invariant 
under rotations in both domain and the target spaces, possessing an 0(3) c j oma j n x 
0(3)target symmetry. As noted earlier, in flat two-dimensional space an additional 
potential term is necessary to ensure the existence of stable finite-size solutions. 
Without it, the repulsive effect of the Skyrme term causes the skyrmions to expand 
indefinitely. In the present model, however, the finite geometry of the sphere acts 
as a stabilizer, so a potential term is not required. Furthermore, stable solutions 
exist even without the Skyrme term. In the latter case, we obtain the well known 
0(3) (or CP 1 ) nonlinear sigma model. 43 

As before, the field <fi in this model is an S 2 — > S 2 mapping, so the relevant 
homotopy group is 7T2(5 2 ) = Z, implying that each field configuration is character- 
ized by an integer topological charge B, the topological degree of the map <fi. In 
spherical coordinates B is given by 

B = ± fdfl W^M , (L15 ) 
4-7T J sin 6 

where dfl = sin^dfld^. 

Static solutions are obtained by minimizing the energy functional 

E -\S + + y f -(^^F) ■ w 

within each topological sector. Before proceeding, it is worthwhile to note that 
setting n = in Eq. (|1.16[) yields the energy functional of the 0(3) nonlinear sigma 
model. The latter has analytic minimal energy solutions within every topological 
sector, given by 

0= (sin f(0)can(B<p), sin f(9)ew.(B<p), cos /((?)), (1.17) 

where f(8) = cos _1 (l— 2(l+(Atan#/2) 2B )~ 1 ) with A being some positive number. 43 
These solutions are not unique, as other solutions with the same energy may be 
obtained by rotating (|1.17|) either in the target or in the domain spaces. The 
energy distributions of these solutions in each sector have total energy Eb = 4ttj3. 

Analytic solutions also exist for the energy functional (|1.16[) with the Skyrme 
term only. They too have the rotationally symmetric form (|1.17|) with f(6) = 9 
and total energy Eb = AttB 2 . They can be shown to be the global minima by the 
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following Cauchy-Schwartz inequality: 



1 f jo<£- ( d s4> x d v 4>)\ 2 ( 1 [ A (I f Ar , / d e <t>xd ip (t> 2 



47T 7 sm6> / V 47r J J V 47r 7 sm # 

(1.18) 

The left-hand side is simply B 2 and the first term in parenthesis on the right-hand 
side integrates to 1. Noting that the second term in the RHS is the Skyrme energy 
(without the k 2 /2 factor), the inequality reads E > 4irB 2 , with an equality for the 
rotationally-symmctric solutions . 

1.3.2. Static solutions 

In general, if both the kinetic and Skyrme terms are present, static solutions of the 
model cannot be obtained analytically. This is with the exception of the B = 1 
skyrmion which has an analytic "hedgehog" solution 

4>[b=i] = (sin cos tp, sin 8 sin ip, cos 0) , (1-19) 

E k 2 
with total energy — = 1 H . 

For skyrmions with higher charge, we find the minimal energy configurations 

by utilizing the full-field relaxation method described earlier. In parallel, we also 

apply the rational map approximation method, originally developed for the 3D 

Skyrme model and compare the results with the relaxation method. Let us briefly 

discuss the rational map approximation method: computing the minimum energy 

configurations using the full nonlinear energy functional is a procedure which is both 

time-consuming and resource-hungry. To circumvent these problems, the rational 

map ansatz scheme has been devised. First introduced by Houghton, Manton and 

Sutcliffe, 44 this scheme has been used in obtaining approximate solutions to the 

3D Skyrme model using rational maps between Riemann spheres. Although this 

representation is not exact, it drastically reduces the number of degrees of freedom 

in the problem, allowing computations to take place in a relatively short amount of 

time. The results in the case of 3D Skyrme model are known to be quite accurate. 

Application of the approximation, begins with expressing points on the base 

q 

sphere by the Riemann sphere coordinate z = tan — e l(p . The complex- valued func- 
tion R(z) is a rational map of degree B between Riemann spheres 

R{z) = ^f\, (1.20) 

where p(z) and q(z) axe polynomials in z, such that max[deg(p), deg(g)] = B, and p 
and q have no common factors. Given such a rational map, the ansatz for the field 
triplet is 

. R + R . R- R 1 - \R\ 2 . , , 

*~ ( l + \R\ 2 ' l + \R\ 2 ' 1 + \R\ 2 >- 1 ' 
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It can be shown that rational maps of degree B correspond to field configurations 
with charge B. 44 Substitution of the ansatz (|I.21j) into the energy functional (|1.16p 
results in the simple expression 

f „ „ K 2 



-E = B + —1 , (1.22) 

4tt 2 ' V 1 



with 



X= — 



I f(l + \z\ 2 



dR 



dz 



(1.23) 



(1 + N 2 ) 5 



4tt7 + 

Note that in the k — -> limit, where our model reduces to the 0(3) nonlinear sigma 
model, the rational maps become exact solutions and the minimal energy value E = 
AitB is attained. Furthermore, the minimal energy is reached independently of the 
specific details of the map (apart from its degree), i.e., all rational maps of a given 
degree are minimal energy configurations in the topological sector corresponding to 
this degree. This is a reflection of the scale- and the rotational invariance of the 
0(3) model. 

In the general case where K ^ 0, the situation is different. Here, minimizing 
the energy (|1.22|) requires finding the rational map which minimizes the functional 
X. As we discuss in the next section, the expression for X given in Eq. (|1.23p is 
encountered in the application of the rational map in the context of 3D skyrmions, 
where the procedure of minimizing X over all rational maps of the various degrees 
has been used. 44-46 Here we redo the calculations, utilizing a relaxation method: to 
obtain the rational map of degree B that minimizes I, we start off with a rational 
map of degree B, with the real and imaginary parts of the coefficients of p(z) and 
q(z) assigned random values from the segment [—1,1]. Solutions are obtained by 
relaxing the map until a minimum of X is reached. 

1.3.3. Relation to the 3D Skyrme model 

In the 3D Skyrme model, the rational map ansatz can be thought of as taken in 
two steps. First, the radial coordinate is separated from the angular coordinates by 
taking the SU(2) Skyrme field U (r, 9, (p) to be of the form 

U(r, 0, <p) = exp(*/(r) <j>{6, <p) ■ tr) , (1.24) 

where er = (ci, 02, 03) are Pauli matrices, f(r) is the radial profile function subject 
to the boundary conditions /(0) = ir and /(oo) = 0, and 4>(8,ip) : S 2 1— > S 2 is a 
normalized vector which carries the angular dependence of the field. In terms of 
the ansatz (jl.24j) . the energy of the Skyrme field is 

E = J 4Tr/' 2 r 2 dr+ / 2(/' 2 + l)sin 2 fdr J ([dg<t>) 2 + _L-(0^) 2 )dfi 

smV dr r (MHf da (1 . 25) 

r 2 J sin 2 8 K ' 
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Note that the energy functional ()1 .25(1 is actually the energy functional of our 
model (|1.16j) once the radial coordinate is integrated out. Thus, our 2D model can 
be thought of as a 3D Skyrme model with a 'frozen' radial coordinate. 

The essence of the rational map approximation is the assumption that <j>(6, ip) 
takes the rational map form (|1.2ip . which in turn leads to a simple expression for 
the energy 



where X is given in Eq. (jl.23p . As in the baby model on the two-sphere, minimizing 
the energy functional requires minimizing X over all maps of degree j3, which is 
then followed by finding the profile function /(r). 

Since the symmetries of the 3D skyrmions are determined solely by the angular 
dependence of the Skyrme field, it should not be too surprising that the solutions 
of the model discussed here share the symmetries of the corresponding solutions of 
the 3D Skyrme model. 

1.3.4. Results 

The configurations obtained from the full-field relaxation method are found to have 
the same symmetries as corresponding multi-skyrmions of the 3D model with the 
same charge. The B = 2 solution is axially symmetric, whereas higher-charge 
solutions were all found to have point symmetries which are subgroups of 0(3). For 
B = 3 and B = 12, the skyrmions have a tetrahedral symmetry. The B = 4 and 
B = 13 skyrmions have a cubic symmetry, and the B = 7 is dodecahedral. The 
other skyrmion solutions have dihedral symmetries. For B = 5 and B = 14 a D^d 
symmetry, for B = 6, 9 and 10 a D^d symmetry, for B = 8 a Ded symmetry and 
for B = 11 a Z?3/j symmetry. In Fig. 11.131 we show the energy distributions of the 
obtained solutions for k 2 = 0.05. 

While for solutions with B < 8 the energy density and the charge density are 
distributed in distinct peaks, for solutions with higher charge they are spread in 
a much more complicated manner. The total energies of the solutions (divided by 
4ttB) are listed in Table along with the symmetries of the solutions (again 

with k 2 = 0.05). 

Application of the rational map ansatz yields results with only slightly higher 
energies, only about 0.3% to 3% above the full-field results. The calculated values 
of X are in agreement with results obtained in the context of 3D skyrmions. For 
9 < B < 14, the rational map approximation yielded slightly less symmetric solu- 
tions than the full-field ones. Considering the relatively small number of degrees of 
freedom, this method all-in-all yields very good approximations. The total energies 
of the solutions obtained with the rational map approximation is also listed in Table 




(1.26) 
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(a)B = 2 (b)B=3 




(c)B=4 (d)B = 5 (e)B=6 




(f)B = 7 (g)B = 8 (h)fl=9 




(1) B = 13 (m) B = 14 



Fig. 1.13. The energy distributions of the multi-skyrmion solutions for charges 2 < B < 14 
(k 2 = 0.05). 



January 12, 2009 1:7 World Scientific Review Volume - 9.75in x 6.5in 



28 



Table 1.1. Total energies (divided by 4irB) of the multi-solitons of the 
model for k 2 = 0.05. 



Charge 


Total energy 


Total energy 


Difference 


Symmetry of 


B 


Full-field 


Rational maps 


in % 


the solution 


2 


1.071 


1.073 


0.177 


Toroidal 


3 


1.105 


1.113 


0.750 


Tetrahedral 


4 


1.125 


1.129 


0.359 


Cubic 


5 


1.168 


1.179 


0.958 


D 2d 


6 


1.194 


1.211 


1.426 


D 4d 


7 


1.209 


1.217 


0.649 


Icosahedral 


8 


1.250 


1.268 


1.406 


D &d 


9 


1.281 


1.304 


1.771 


D 4d 


10 


1.306 


1.332 


1.991 


D 4d 


11 


1.337 


1.366 


2.224 


D ih 


12 


1.360 


1.388 


2.072 


Tetrahedral 


13 


1.386 


1.415 


2.137 


Cubic 


11 


1.421 


1.459 


2.712 


D 2 



1.3.5. Further remarks 

As we have just seen, the baby Skyrme model on the two-sphere shares very signif- 
icant similarities with the 3D model, especially in terms of multi-skyrmion symme- 
tries. The fact that the model discussed here is two-dimensional makes it simpler to 
study and perform computations with, when compared with the 3D Skyrme model. 

Some of the results presented above may, at least to some extent, also be linked to 
the baby skyrmions which appear in two-dimensional electron gas systems, exhibit- 
ing the quantum Hall effect. As briefly noted in the Introduction, baby skyrmions 
arise in quantum Hall systems as low-energy excitations of the ground state, near fer- 
romagnetic filling factors (notably 1 and 1/3). 6 ' 47 There, the skyrmion is a twisted 
two-dimensional configuration of spin, and its topological charge corresponds to the 
number of time the spin rotates by 2tt. While for the electron gas, the stability 
of the soliton arises from a balance between the electron-electron Coulomb energy 
and the Zeeman energy, in our model the repulsive Skyrme-term energy is balanced 
by the underlying geometry (i.e., the sphere). The connection between these two 
models suggests possible existence of very structured spin textures in quantum Hall 
systems, although a more detailed analysis of this analogy is in order. 

1.4. Rotating Baby Skyrmions 

We now turn to analyze the phenomenon of spontaneous breaking of rotational 
symmetry (SBRS) as it appears in rotating baby skyrmions. In general, SBRS refers 
to cases where physical systems which rotate fast enough deform in a manner which 
breaks their rotational symmetry a symmetry they posses when static or rotating 
slowly. The recognition that rotating physical systems can yield solutions with less 
symmetry than the governing equations is not new. One famous example which 
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dates back to 1834 is that of the equilibrium configurations of a rotating fluid mass. 
It was Jacobi who was first to discover that if rotated fast enough, a self-gravitating 
fluid mass can have equilibrium configurations lacking rotational symmetry. In 
modern terminology Jacobi's asymmetric equilibria appear through a symmetry 
breaking bifurcation from a family of symmetric equilibria as the angular momentum 
of the system increases above a critical value (a "bifurcation point"). 48 ' 49 Above 
this critical value, rotationally-symmetric configurations are no longer stable, and 
configurations with a broken rotational symmetry become energetically favorable. 

By now it is widely recognized that symmetry-breaking bifurcations in rotating 
systems are of frequent occurrence and that this is in fact a very general phe- 
nomenon, appearing in a variety of physical settings, among which are fluid dy- 
namics, star formation, heavy nuclei, chemical reactions, plasmas, and biological 
systems, to mention some diverse examples. 

Recently, SBRS has also been observed in self-gravitating iV-body systems, 50 ' 51 
where the equilibrium configurations of an iV-body self-gravitating system enclosed 
in a finite 3 dimensional spherical volume have been investigated using a mean- 
field approach. It was shown that when the ratio of the angular momentum of the 
system to its energy is high, spontaneous breaking of rotational symmetry occurs, 
manifesting itself in the formation of double-cluster structures. These results have 
also been confirmed with direct numerical simulations. 52 

It is well-known that a large number of phenomena exhibited by many-body 
systems have their counterparts and parallels in field theory, which in some sense is 
a limiting case of iV-body systems in the limit N — > oo. Since the closest analogues 
of a lump of matter in field theories are solitons, the presence of SBRS in self- 
gravitating TV-body systems has led us to expect that it may also be present in 
solitonic field theories. 

Our main motivation towards studying SBRS in solitons is that in hadronic 
physics Skyrme-type solitons often provide a fairly good qualitative description 
of nucleon properties (See, e.g., 53,54 ). In particular, it is interesting to ask what 
happens when such solitons rotate quickly, because this might shed some light on the 
non-spherical deformation of excited nucleons with high orbital angular momentum, 
a subject which is now of considerable interest. 

We shall see that the baby Skyrme model on the two-sphere indeed exhibits 
SBRS, and we will try to understand why this is so. 55 First, we give a brief account 
for the occurrence of SBRS in physical systems in general, and then use the insights 
gained from this discussion to infer the conditions under which SBRS might appear 
in solitonic models and in that context we study its appearance in baby Skyrme 
models. Specifically, we shall show that SBRS emerges if the domain manifold of 
the model is a two-sphere, while if the domain is R 2 , SBRS does not occur. 
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1.4.1. SBRS from a dynamical point of view 

The onset of SBRS may be qualitatively understood as resulting from a competition 
between the static energy of a system and its moment of inertia. To see this, let us 
consider a system described by a set of degrees of freedom cj), and assume that the 
dynamics of the system is governed by a Lagrangian which is invariant under spatial 
rotations. When the system is static, its equilibrium configuration is obtained by 
minimizing its static energy -E s tatic with respect to its degrees of freedom <f> 

SE 

— = where E = £ static (0) . (1.27) 

Usually, if £ , static(</') does not include terms which manifestly break rotational sym- 
metry, the solution to (|1.27|) is rotationally-symmetric (with the exception of de- 
generate spontaneously-broken vacua, which are not of our concern here). If the 
system rotates with a given angular momentum J = Jz, its configuration is natu- 
rally deformed. Assuming that the Lagrangian of the system is quadratic in the time 
derivatives, stable rotating configurations (if such exist) are obtained by minimizing 
its total energy Ej 

8E J 2 

— f = where £ 7 = £ static (0) + — — , (1.28) 
d<p 2I((j>) 

where /(</>) is the ratio between the angular momentum of the system and its angular 
velocity u) = ujz (which for simplicity we assume is oriented in the direction of the 
angular momentum). I ((f) is the (scalar) moment of inertia of the system. 

The energy functional (|1.28[) consists of two terms. The first, ©static, increases 
with the asymmetry. This is simply a manifestation of the minimal-energy config- 
uration in the static case being rotationally-symmetric. The second term J 2 /2I, 
having the moment of inertia in the denominator, decreases with the asymmetry. 
At low values of angular momentum, the -©static term dominates, and thus asymme- 
try is not energetically favorable, but as the value of angular momentum increases, 
the second term becomes dominant, giving rise to a possible breaking of rotational 
symmetry. 



1.4.2. SBRS in baby Skyrme models 

In what follows, we show that the above mechanism of SBRS is present in certain 
types of baby Skyrme models. 

As already discussed in previous sections, the static solutions of the baby Skyrme 
model have rotationally-symmetric energy and charge distributions in the 

charge-one and charge-two sectors. 4 The charge-one skyrmion has an energy peak 
at its center which drops down exponentially. The energy distribution of the charge- 
two skyrmion has a ring-like peak around its center at some characteristic distance. 
The rotating solutions of the model in R 2 are also known. 15,56 Rotation at low angu- 
lar velocities slightly deforms the skyrmion but it remains rotationally-symmetric. 
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For larger values of angular velocity, the rotationally-symmetric configuration be- 
comes unstable but in this case the skyrmion does not undergo symmetry breaking. 
Its stability is restored through a different mechanism, namely that of radiation. 
The skyrmion radiates out the excessive energy and angular momentum, and as a 
result begins slowing down until it reaches equilibrium at some constant angular 
velocity, its core remaining rotationally-symmetric. Moreover, if the Skyrme fields 
are restricted to a rotationally-symmetric (hedgehog) form, the critical angular ve- 
locity above which the skyrmion radiates can be obtained analytically. It is simply 
the coefficient of the potential term u cr it = /J.. 15 Numerical full-field simulations 
also show that the skyrmion actually begins radiating well below cj cr it, as radiation 
itself may be non-rotationally-symmetric. The skyrmion's core, however, remains 
rotationally-symmetric for every angular velocity. 

The stabilizing effect of the radiation on the solutions of the model has lead 
us to believe that models in which radiation is somehow inhibited may turn out 
to be good candidates for the occurrence of SBRS. In what follows, we study the 
baby Skyrme model on the two-sphere, whose static solutions were presented in 
the previous section. Within this model, energy and angular momentum are not 
allowed to escape to infinity through radiation, and as a consequence, for high 
enough angular momentum the mechanism responsible for SBRS discussed in the 
previous section takes over, revealing solutions with spontaneously broken rotational 
symmetry. 

1.4.3. The baby Skyrme model on the two-sphere 

In order to find the stable rotating solutions of the model, we assume for simplicity 
that any stable solution would rotate around the axis of angular momentum (which 
is taken to be the z direction) with some angular velocity uj. The rotating solutions 
thus take the form </>(#, ip, t) = <f>{6, if — ut). The energy functional to be minimized 
is 



where / is the ratio of the angular momentum of the skyrmion to its angular velocity, 
or its "moment of inertia" , given by 



1.4.4. Results 

In what follows we present the results obtained by the minimization scheme applied 
to the rotating solutions of the model in the charge-one and charge-two sectors, 
which as mentioned above are rotationally-symmetric. For simplicity, we fix the 
parameter k at k 2 = 0.01 although other k values were tested as well, yielding 
qualitatively similar solutions. 



E — -Ecstatic + 777 i 



(1.29) 




(1.30) 
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I. 4.4.1. Rotating charge-one solutions 

The rotating charge-one skyrmion has spherically-symmetric energy and charge dis- 
tributions in the static limit (Fig. 11.14b ,). When rotated slowly, its symmetry is 
reduced to 0(2), with the axis of symmetry coinciding with the axis of rotation (Fig. 

II. 14b ) . At some critical value of angular momentum (which in the current settings 
is Jcrit ~ 0.2), the axial symmetry is further broken, yielding an ellipsoidal energy 
distribution with three unequal axes (Fig. 11.14b ). Any further increase in angular 
momentum results in the elongation of the skyrmion in one horizontal direction and 
its shortening in the perpendicular one. The results arc very similar to those of the 
rotating self-gravitating ellipsoid. 




(a) J = (b) J = 0.14 (c) J = 0.26 

Fig. 1.14. Baby skyrmions on the two-sphere (ft 2 = 0.01): The charge distribution B(9, ip) of the 
charge-one skyrmion for different angular momenta. In the figure, the vector B(8, tp)r is plotted 
for the various 6 and tp values. 
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1.4.4.2. Rotating charge-two solutions 

SBRS is also observed in rotating charge-two skyrmions. The static charge-two 
skyrmion has only axial symmetry (Fig. 11.15a ). with its symmetry axis having 
no preferred direction. Nonzero angular momentum aligns the axis of symmetry 
with the axis of rotation. For small values of angular momentum, the skyrmion is 
slightly deformed but remains axially-symmetric (Fig. 11.15b ). Above J cr it ~ 0.55 
however, its rotational symmetry is broken, and it starts splitting to its 'constituent' 
charge-one skyrmions (Fig. 11.15b and ll,15b l). As the angular momentum is further 
increased, the splitting becomes more evident, and the skyrmion assumes a string- 
like shape. This is somewhat reminiscent of the well-known elongation, familiar 
from high-spin hadrons which are also known to assume a string-like shape with the 
constituent quarks taking position at the ends of the string. 58,59 




(a) J = o 



(b) J = 0.55 



(c) J = 0.1 





(d) J = 0.66 



(e) J = 0.75 



Fig. 1.15. Baby skyrmions on the two-sphere (ft 2 = 0.01): The charge distribution B(9, ip) of the 
charge-two skyrmion for different angular momenta. In the figure, the vector B(8, ip)v is plotted 
for the various 6 and ip values. 



1.4.5. The rational map ansatz 

A somewhat more analytical analysis of this system may be achieved by the use 
of the rational maps approximation scheme, 44 which as was shown earlier provides 
quite accurate results for the static solutions of the model. 57 

In its implementation here, we simplify matters even more and reduce the de- 
grees of freedom of the maps by a restriction only to those maps which exhibit the 
symmetries observed in the rotating full-field solutions. This allows the isolation 
of those parameters which are the most critical for the minimization of the energy 
functional. 
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As shown in Fig. 11.141 the charge and energy densities of the charge-one skyrmion 
exhibit progressively lower symmetries as J is increased. The static solution has an 
0(3) symmetry, while the slowly-rotating solution has an 0(2) symmetry. Above a 
certain critical J, the 0{2) symmetry is further broken and only an ellipsoidal sym- 
metry survives. Rational maps of degree one, however, cannot produce charge den- 
sities which have all the discrete symmetries of an ellipsoid with three unequal axes. 
Nonetheless, approximate solutions with only a reflection symmetry through the xy 
plane (the plane perpendicular to the axis of rotation) and a reflection through one 
horizontal axis may be generated by the following one-parametric family of rational 
maps 

, , cos a 
R{z) = ; , 1.31 

which has the charge density 

B(9,<p)=( r^~ Q V- (!- 32 ) 

\ 1 + sin a sin cos ip J 

Here, a 6 [— 7r, tt] is the parameter of the map, with a = corresponding to 
a spherically-symmetric solution and a non-zero value of a corresponding to a 
nonrotationally-symmetric solution. Results of a numerical minimization of the 
energy functional (| 1 . 2Q[) for fields constructed from (]1.3ip for different values of 
angular momentum J are shown in Fig. 11.16b . While for angular momentum less 
than J cr ;t k 0.1, a = minimizes the energy functional (a spherically-symmetric 
solution), above this critical value bifurcation occurs and a = is no longer a mini- 
mum; the rotational symmetry of the charge-one skyrmion is broken and it becomes 
nonrotationally-symmetric. 

A similar analysis of the charge-two rotating solution yields the one-parametric 
map 

. sin a + z 2 cos a .„ „„. 

R(z) = , 1.33 

cos a + z £ sin a 

with corresponding charge density 

»<'•«" -frr^T^V-n)'- (L34) 

\2 + sin t/(sm 2a cos 2ip — 1) / 

In this case, a = corresponds to a torodial configuration, and a non-zero value of a 
yields solutions very similar to those shown in Fig. I1.15[ having the proper discrete 
symmetries. The results in this case are summarized in Fig. 11.16b . indicating that 
above J cr it ~ 0.57 the minimal energy configuration is no longer axially-symmetric. 

The discrepancies in the critical angular momenta J cr ;t between the full-field 
method (0.2 for charge-one and 0.55 for charge-two) and the rational maps scheme 
(0.1 for charge-one and 0.57 for charge-two) are of course expected, as in the latter 
method, the solutions have only one degree of freedom. Nonetheless, the qualitative 
similarity in the behavior of the solutions in both cases is strong. 



January 12, 2009 1:7 World Scientific Review Volume - 9.75in x 6.5in text 



30 M. Karliner and I. Hen 



0.05 0.10 0.15 0.3 0.6 




0.05 0.10 0.15 0.3 0.6 

J J 

(a) (b) 



Fig. 1.16. Spontaneous breaking of rotational symmetry in the restricted rational maps approx- 
imation for the baby skyrmions on the two-sphere: the parameter a as a function of the angular 
momentum J, for the charge-one (top) and the charge-two (bottom) solutions. The lines are to 
guide the eye. 

1.4.6. Further remarks 

We have seen that SBRS appears not only in rotating classical-mechanical systems 
but also in the baby Skyrmc model on the two-sphere. We have argued that this 
is so because the phenomenon originates from general principles, and hence it is a 
universal one. 

The results presented above may, at least to some extent, also be linked to recent 
advances in the understanding the non-sphericity of excited nucleons with of large 
orbital momentum. Non-spherical deformation of the nucleon shape is now a focus 
of considerable interest, both experimental 60 ' 61 and theoretical. 62-64 As skyrmions 
are known to provide a good qualitative description of many nucleon properties, 
the results presented here may provide some corroboration to recent results on this 
subject (e.g., 64 ), although a more detailed analysis of this analogy is in order. 
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A.l. Obtaining Baby Skyrmion Solutions — The Relaxation 
Method 

As a large part of the studies presented above is based on numerically obtaining the 
baby skyrmion configurations, in the following we describe the relaxation method 
that was used to obtain the solutions. 

The multi-solitons of the baby Skyrme model are those field configurations which 
minimize the static energy functional within each topological sector. The energy 
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functional is given by 

E= jd 2 X Q(9^ • d X <f> + dy<P ■ dy + y X dyCj)) 2 + Ufa)) . (A.l) 

As already noted, the baby Skyrme model is a nonintegrable system, so in 
general, explicit analytical solutions to its Euler-Lagrange equations are nearly im- 
possible to find. Hence, one must resort to numerical techniques. 

Generally speaking, there are two main approaches to finding the baby skyrmion 
solutions numerically. One approach is to employ standard techniques to numer- 
ically solve the Euler-Lagrange equations which follow from the energy functional 
(|A.1|) . The other approach - the one taken here - is to utilize relaxation methods 
to minimize the energy of the skyrmion within any desired topological sector. In 
what follows, we describe in some detail the relaxation method we have used all 
throughout this research. This method is based on the work of Hale, Schwindt and 
Weidig. 67 We assume for simplicity that the base space is descretized to a rectan- 
gular grid. The implementation of this method in the case of curved spaces or for 
a polar grid is straightforward. 

The relaxation method begins by defining a grid with N 2 points, where at each 
point a field triplet 4>(x mi y n ) is defined. All measurable quantities such as energy 
density or charge density are calculated at the centers of the grid squares, using the 
following expressions for the numerical derivatives, also evaluated at these points: 



90 

dx 



= J- x (A.2) 

O m+ i,y„+i) Ax 



/ / 4>(x m+1 , y n ) + 4>{x m +l,Vn+l) \ / 4>{Xm,Vn) + <j>{x m , Vn+l) ' 

\ \ 2 / normcd \ 2 / normed 

with the y-derivatives analogously defined, and the "normed" subscript indicates 
that the averaged fields are normalized to one. If the field itself has to calculated 
at that center of a grid square, we use the prescription 

(p(x m+ i,y n+ i) (A.3) 

= (- (cf)(x m ,y n ) + <p(x m ,y n+1 ) + 4>{x m +\,y n ) + cf)(x m+1 ,y n+1 )) ) 

\ 4 / normed 

The basic updating mechanism of the relaxation process consists of the following 
two steps: A point (x m ,y n ) on the grid is chosen at random, along with one of the 
three components of the field (p(x m ,y n ). The chosen component is then shifted by 
a value chosen uniformly from the segment [— A^, A^] where A^ = 0.1 initially. 
The field triplet is then scaled and the change in energy is calculated. If the energy 
decreases, the modification of the field is accepted and otherwise it is discarded. 

The relaxation process, through which the energy of the baby skyrmion is min- 
imized, is as follows: 
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(1) 



Initialize the field triplet <fi to a rotationally-symmetric configuration 



^initial = (sin f(r) cos B8, sin f(r) sin B6, cos /(r)) 



(A.4) 



In our setup, we have chosen the profile function /(r) to be f(r) = Trexp(— r), 
r and 9 being the usual polar coordinates. 

(2) Perform the basic updating mechanism for M x TV 2 times (we took M = 100), 
and then calculate the average rate of acceptance. If it is smaller than 5%, 
decrease by half. 

(3) Repeat step ^ until < 10~ 9 , meaning no further decrease in energy is 
observed. 

This procedure was found to work very well in practice, and its accuracy and 
validity were verified by comparison of our results to known ones. There is however 
one undesired feature to this minimization scheme, which we note here: it can get 
stuck at a local minimum. This problem can be resolved by using the "simulated 
annealing" algorithm, 65,66 which in fact has been successfully implemented before, 
in obtaining the minimal energy configurations of three dimensional skyrmions. 67 
The algorithm is comprised of repeated applications of a Metropolis algorithm with 
a gradually decreasing temperature, based on the fact that when a physical system 
is slowly cooled down, reaching thermal equilibrium at each temperature, it will 
end up in its ground state. This algorithm, however, is much more expensive in 
terms of computer time. We therefore employed it only in part, just as a check 
on our method, which corresponds to a Metropolis algorithm algorithm at zero 
temperature. We found no apparent changes in the results. 
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